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Abstract 


The present research study describes the computational P-matrix model based modelling and 
subsequent analysis of the frequency response of a 300 MHz ST-X Quartz Surface Acoustic Wave (SAW) 
delay line device designed and fabricated with 43.5 IDT finger pairs. Employing a custom made 
MATLAB® algorithm, the device is modelled with a provision to vary its design parameters. Finally, a 
comparison of the computationally simulated frequency response with experimental device results 
indicates promising and good agreement between model and experiment which can be utilized and 
implemented during design and development stage of these devices for their potential applications in 


SAW Sensor systems. 
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1. INTRODUCTION 


Surface Acoustic Wave (SAW) sensors (Raj 
et al. 2013) are extremely versatile devices that are 
commercially possible for transduction of physical 
and chemical quantities. Almost all physical quantities 
like Surface mass (HribSek et al. 2010), Stress 
(Filipiak et al. 2007), Strain (Hempel ert al. 2013), 
temperature(Malocha et al. 2013), Pressure (Benetti et 
al. 2008) etc., can be measured by SAW sensors. 
SAW sensors are miniaturized, highly sensitive, 
reproducible, long term stable and fast real time 
response (Banu Priya et al. 2014). 


Computer models for SAW devices are 
continuously improved to meet the demanding design 
requirements for key components in _— sensor 
applications(Amani et al. 2012; Chu et al. 2013; 
Sharma et al. 2014). This often leads to the increase in 
simplicity of the simulation codes and reduce in run 
time of the program. Recently, some sophisticated 
models have been developed. In our present study, a 
straight forward approach like P-Matrix model 1s 
employed (Rusakov et al. 2001). It provides fast and 
accurate simulation, which are used for designing 
SAW devices (Pastureaud, 2004). 


Modelling of SAW devices is resorted to 
achieve two aims : 
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(1) to understand propagation, generation and 
detection of surface waves in piezoelectric materials 
and 

(11) to study and design the structures such 
as IDTs, delay lines, filters, resonators, sensors etc., 
to achieve preferred device responses. 


Modelling techniques need to provide fast 
and exact results, be comfortable to compute and 
distinctly offer relation between modelled parameters 
and device performance. A ST-X Quartz based 300 
MHz SAW delay line with uniform IDTs possessing 
43.5 finger pairs per IDT has been generated by using 
a MATLAB® software algorithm. Modelled device 
parameters are assumed, studied and _ graphically 
represented. The whole purpose of this simulation 
study is to model the SAW delay line for enhanced 
device performance and subsequently compare the 
results obtained from the experimental one (Haresh et 
al. 2010). 


2. THEORETICAL ANALYSIS 
2.1 P-Matrix modelling 

The P-Matrix has been invented by Tobolka 
in 1979 (Manner et al. 1996) and is now widely used 


(Ruppel et al. 1994; Brown et al. 1989; Ventura et al. 
1992) (Hartmann and Abbott, 1988). It is a mixed 
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matrix representation of an IDT with three ports 
(Fig. 1). This three port has one electrical and two 
acoustical ports. To describe a whole transducer, the 
standard cell (Smith et al. 1972) three ports have to be 
cascaded acoustically whereas the electrical ports 
must be connected in parallel. The standard cell 
originally proposed by Smith et al. 1972 (Tancrell et 
al. 1969; Tancrell and Holland, 1971) represents one 
single finger and is described by means of the crossed- 
field equivalent circuit. Therefore, transducers are 
represented here by a mixed matrix, so called P- 
matrix, defined in this correspondence. The voltage- 
to-SAW transfer functions, the electrical admittance 
and the SAW scattering matrix of the shorted 
transducer are elements of the P-matrix, which will be 
calculated for a standard cell and the whole 
transducer. The two acoustic ports are described by an 
S-Matrix, and the electric port is described by an 
admittance. 
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Fig. 1: Diagram showing mixed matrix variables at each 
port 


P-matrix model is a distinct form of COM 
model, a_ well-established tool to consider the 
electrical and acoustical properties of IDT’s and 
reflector gratings (Tobolka, 1979). This matrix is 
convenient for analysing a wide variety of SAW 
devices. This model relates both outgoing acoustic 
surface waves A,,/A, and the electric current I to both 
incoming acoustic surface waves Aj/Aj. and the 
electric voltage V; (shown in Fig. 2). 





Fig. 2: P-Matrix Parameters of transducer 


The P-matrix is defined such that, 


Pi, Pi2 P13 Ail 
Aj? 


Pz; Px Pr3 2 
P3; P32 P33 Vi 


P - matrix 
(1) 


Here, P33 denotes the transducer admittance 
Y, (scalar), Py7, Pi2, P2; and P22 are scattering 
parameters (matrices), P;; and P23 are electro-acoustic 
excitation matrix (vectors) and P3, and P3 are 
acousto-electric detection matrix (vectors). The value 
of Pj. and P>, elements are non-zero from the basic 
cell idea introduced by (Ruppel et al. 1994). The most 
complicated one is the excitation matrix, where all P;;, 
elements are non-zero. The parameters of the above P- 
Matrix are given below 


3 Ay, = jp,(k)(@W|s/2 exp(—jk L/2) 
P,= Pf, = jp, blow s/2exp(—jk 1/2) 


P., =Y,(@)=G,(@)+ JB, (@)+ ja, 


ao 
lI 


(2) 


From the above equations, we define a 
function —, (x) as the electrostatic charge density on 


the electrodes when unit voltage is applied and is 
calculated by a Green's function method (Baghai- 
Wadji et al. 1984). kis the wave number for waves 
between the two port location, L is the length of the 


IDT. The constant Is is defined as 


 _ (Av/¥) (2) 
é 


oe) 


From the above equation €, is known as 
effective permittivity (Guerra-Pulido and _ Pérez- 
Alcazar, 2014). The conductance G,(@) can be 


written as 
=. 2 
G,(@) = oWls|p,(b)} 3) 
Where, @ is angular frequency and W is the 


aperture of uniform IDT. The susceptance B, (@) can 


be calculated from the equation, 
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The insertion loss for the SAW device can be 
determined by the quantity -20 log J, in dB,is 
2G.G 
=F, i LS —_ 2 
IPs +Y,, 


(5) 


(6) 
Where G, = Re 1Y ae P, and P; are source 


power and power deliver to the load. 
2.1 Cascading of two port components 


In the present study, we considered two 
adjacent transducers with uniform aperture W and P- 


for SAW devices as shown in 


matrices P and Ee 
Fig. 3. 





Fig. 3: P-Matrix for two component device 
For a symmetric device structure, the 


admittance matrix elements (Y-parameters) of the two 
port SAW delay line can be written as 


2p) (Pi) ei 


Y, =Pi- 

a 
-2(P3P5)e™ 

Ys Ve 

12 21 1 (PcP? Je (3) 


Where Y;; and Y2; represent the input and 
output (transfer) admittance(Guerra-Pulido and Pérez- 
Alcazar, 2014). Using the admittance matrices, the 
transfer function of the device can be obtained. For a 
two-port resonator device is connected to an arbitrary 
source and load, then by using Y to S conversion 
relations, the device transfer function can be written as 


2Y,,VR,R, 


Re $Y ROG Yok) 


Mee eo e, 


12 


S,=S 


(9) 
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The above equation gives the complete 
response of the SAW devices that has been derived by 
by P-matrix. 


2.2 Modelling strategy 


The various steps adapted for modelling of SAW 
devices are highlighted as follows: 

1. In present study ST-X Quartz substrate was 
selected. Because compare’ to __ other 
substrates,quartz has a good _ temperature 
stability,the temperaturs coefficient of delay 
(TCD) is nearly zero(Campbell, 1989). 

IDT material chosen is Aluminium (A\). 
Because Al serves as an inert metal as well as 
good adhesive( Venkatesan and Pandya, 2013). 
The input parameters for SAW _ device 
modelling according to P-matrix model is 
shown in Table 1. 


Table 1. Input parameters for modelling of SAW Device 





S. No Parameter (Symbol) Values 


2 0.0016 (ST-X 


| 2 | saw Velocity (Vs) 3158 m/s 


= Operating frequency (f)) | 300 MHz 


IDT geometry Split geometry 
s Number of finger pairs 1a aeons 

(Np) 

Load and Source 
Employing the P- Matrix model, the outputs are 
shown in graphical format to compare 
experimental results. 
An operating center frequency of 300 MHz was 
modelled using a MATLAB algorithm and 
simulated results of SAW device output graphs 


were obtained. The modelling strategy adapted 
as showing Fig. 7. 





Input Parameters 


MATLAB Program 


Output Values/Graphs_ 


Comparative Result 





Fig. 7: Flow chart of Modelling Strategy 
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3. RESULT & DISCUSSION 


As previously declared, an ST-X Quartz 
piezoelectric substrate based 300MHz SAW delay line 
with 43.5 finger pairs, split electrode geometry and 
constant aperture was fabricated in-house employing 
laser pattern writing for mask making, thermal 
evaporation for metallization and Photolithographic 
(PLG) etching method for device making. This was 
followed by the standard procedures of dicing, 
bonding and packaging in the form of a portable unit 
(Haresh M. Pandya et al. 2013). 


The modelled device response was calculated 
as given by equation (9). The steps involved in the 
modelling procedure were codified into an intensely 
computative MATLAB® program which was run 
many times by varying the free surface SAW velocity 
(vo) for ST-X Quartz piezoelectric substrate around 
3158 m/s. The results are also exported to .x/s format 
for comparison with the experimental values obtained. 
The MATLAB® algorithm was codified to further 
generate plots of modelled parameters as a function of 
frequency so as to deduce their behavior directly. 


In addition to the above, the algorithm also 
had the flexibility to simulate device parameters for 
various sets of SAW delay lines for the device 
designer. The striking point of this modelling exercise 
was that the specific frequency response was obtained 
for the 300 MHz SAW delay line. The modelled 
parameters received from the P-Matrix Model are 
summarized in Table 2. 


Measured Response 
—— Modelled Response 





-29.2029 dB 
-29.5592 dB 


Hil 


310.00 


i 


280.00 


285.00 290.00 295.00 300.00 315.00 320.00 


Fig. 8. Comparison of Modelled and Measured Values 


A comparison of the modelled frequency 
response with the measured values (Fig. 8) 
demonstrates in an excellent manner, the near perfect 
match achieved in the frequency response plots which 
account for second order effects such as EMF, TTE 
and circuit parasitics. All important parameters like 
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observed frequency response, insertion loss and 3dB 
bandwidth are accurately modelled and predicted. A 
difference of around 0.3583 dB is observed between 
the modelled value of insertion loss -29.5592 dB and 
its measured value of -29.2029 dB, the difference of 
around 0.3660 MHz is observed between the modelled 
value of 3dB bandwidth 4.1315 MHz and its measured 
value of 4.4975 MHz and frequency responses also 
nearly same. 


It canbe clearly seen that there is a slight 
difference in the modelled values of IL as shown in 


Table 2. Comparison of Modelled and Measured Values 


S. No Parameter Modelled Measured 
(symbol) Values Values 


Operating 299.99 299.59 
frequency (fo) MHz MHz 
o oe -29.2029 
Insertion Loss 
dB 
ETE 4.4975 
3 3dB Bandwidth MHz MHz 





4. CONCLUSION 


The present research study describes the 
computational P-matrix model based modelling and 
subsequent analysis of the frequency response of a 
300 MHz ST-X Quartz Surface Acoustic Wave 
(SAW) delay line device designed and fabricated with 
43.5 IDT finger pairs. 


Employing a custom made MATLAB® 
algorithm, the device is modelled with a provision to 
vary its design parameters. Finally, a comparison of 
the computationally simulated frequency response 
with experimental device results indicates promising 
and good agreement between model and experiment 
which can be utilized and implemented during design 
and development stage of these devices for their 
potential applications in SAW Sensor systems. 
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